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Abstract 

A quasipotential formalism for elastic scattering from relativistic bound states 
is based on applying an instant constraint to both initial and final states in 
the Breit frame. This formalism is advantageous for the analysis of elec- 
tromagnetic interactions because current conservation and four momentum 
conservation are realized within a three-dimensional formalism. Q Wave func- 
tions are required in a frame where the total momentum is nonzero, which 
means that the usual partial wave analysis is inapplicable. In this work, the 
three-dimensional equation is solved numerically, taking into account the rel- 
evant symmetries. A dynamical boost of the interaction also is needed for the 
instant formalism, which in general requires that the boosted interaction be 
defined as the solution of a four-dimensional equation. For the case of a scalar 
separable interaction, this equation is solved and the Lorentz invariance of the 
three-dimensional formulation using the boosted interaction is verified. For 
more realistic interactions, a simple approximation is used to characterize the 
boost of the interaction. 

25.30.Bf, 24.10.Jv, ll.10.St 
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I. INTRODUCTION 



The theory of relativistic bound states in quantum field theory features four- dimensional 
equations and, in principle, an infinite number of degrees of freedom. Practical methods 
to solve the problem are not available and it is therefore common to reduce its complexity. 
In this paper, we discuss a covariant reduction to three dimensions and a finite number of 
degrees of freedom. |lj The interactions used are instantaneous. The formalism is tractable 
and it is applied to the deuteron bound state problem. 

A general technique for reducing four- dimensional dynamics to three dimensions is to 
introduce a constraint which fixes one component of a four-vector in terms of the others. ]2[ 
This may be done covariantly, but the resulting formalism may possess unphysical asymme- 
tries or singularities. An example is when one particle is constrained to its mass shell, and a 
second, identical particle is not. |J Symmetry with respect to exchange of particle labels is 
lost. It is possible to respect the Pauli Principle by use of appropriately defined interactions, 
but they are cumbersome and they possess unphysical singularities which must be removed 
by hand. f| 

Symmetrical three-dimensional reductions have been used recently by Tjon and collabo- 
rators. |5||| Our work is similar but is based on the three-dimensional formalism developed 
by Mandelzweig and Wallace , with an instant constraint at non-zero total momentum. 
IHII An instant constraint maintains symmetry with respect to exchange of particle labels and 
yields a constrained equation for relativistic bound states with no unphysical singularities. 
In this regard, it is an attractive alternative to reductions in which one particle is constrained 
to its mass shell. Although unphysical singularities are avoided, they are not entirely absent. 
They are expected to appear in corrections to the theory and in the four-dimensional equa- 
tions for the boost of the interaction. Unphysical singularities are an unsolved problem for 
quasipotential approaches. In our approach to elastic electron scattering from a two-body 
bound state, they play no role. 

A theory of bound states needs to be complemented by a corresponding theory of cur- 
rents, e.g., the electromagnetic current. It is necessary to maintain a well-defined connection 
to quantum field theory and, for this reason, the Bethe-Salpeter formalism is the preferred 
starting point for the two-body problem. The associated currents have been formulated by 
Mandelstam in a celebrated paper. In general, the currents depend on the interactions 
used to describe the bound state and a consistent relationship between the two is manda- 
tory in order to conserve the electromagnetic current. An attractive feature of the instant 
formalism we discuss is that a Ward-Takahashi identity is realized at the level of the rel- 
ativistic impulse approximation, thus guaranteeing current conservation for elastic matrix 
elements. There is a price associated with the current conservation. In the analysis of elastic 
electromagnetic scattering, a boost of the interaction is necessary in order to calculate the 
bound state wave functions in the Breit frame, where the total three-momentum is nonzero. 

Dirac showed that the boost operator must depend on the interactions when the genera- 
tors of the Poincare group are quantized at an instant of time. |TU| This has its counterpart 



in the instant quasipotential formalism. A dynamical, four-dimensional equation must be 
solved to determine the instant quasipotential corresponding to different values of the total 
three-momentum. We distinguish between a kinematical boost and a dynamical one. The 
former is simpler because it is effected in the same fashion as for free particles and this is 
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an advantage of a covariant formalism. [§,|J. The latter is required in an instant formalism. 
In the work of Hummel and Tjon, || an instant quasipotential is used together with a kine- 
matical boost. This leads to an inconsistency with respect to current conservation in matrix 
elements. In the present work, we develop methods to handle the dynamical boost that are 
consistent with current conservation for electromagnetic interactions. 

For the special case of a scalar, separable interaction, the dynamical equation correspond- 
ing to a boost of the interaction is solved exactly and we verify that the mass of the deuteron 
is invariant when it is used. The boosted interaction varies rather slowly with the momen- 
tum for the deuteron and it may be approximated by a renormalization of the strength of 
the rest-frame interaction. This approximation is used in our analysis of electromagnetic 
form factors, which will be the subject of another paper. 

This work contains the analysis and methods of solution for wave functions in the Breit 
frame. Because the total three-momentum is nonzero, the usual partial wave analysis is 
inapplicable. The wave functions are obtained by solving the three-dimensional integral 
equations using appropriately boosted interactions. Section II reviews the quasipotential 
equations in a frame where total three-momentum is nonzero. Section III discusses the 
equation which relates the quasipotential in different frames and presents the solution for 
a scalar, separable interaction. In Section IV, we discuss the symmetries of the relativistic 
bound state equations and use them to reduce the equations to a solvable form for the case 
of a boson-exchange interaction. Section V presents results for the deuteron wave functions 
and their variation with the total momentum of the bound state. Section VI presents some 
concluding remarks. 



II. RELATIVISTIC BOUND-STATE EQUATION 

The two-body equation with instant constraint has been derived in Ref. in the rest 
frame of the two-body system. The derivation incorporates crossed graphs using a form 
of the eikonal approximation. The same derivation carried out in a frame where the total 
momentum is nonzero yields the three-dimensional quasipotential equation, 

—V(p,k; P)m;P), (1) 

with relative and total momenta p = {jp-y — £>2)/2 and P = p\ + P2 (p = 0, P = P z z). The 
three-dimensional propagator is expressed in terms of projection operators for positive- and 
negative-energy states as follows, 

a (o P) - V Af(Pi)A?(p 2 ) m 
Sl(P ' P) -^(p 1+ P 2 )(Bp/2)- ei - £2 ' (2) 



where pi = ±, Af (p$) = piPf 1 ^, where Pf* = wf (pi'Pi)[Ui i (piPi)y are projection operators 
obeying P^P?* = 5 Pi:P '_. Dirac spinors used obey the hermitian normalization condition (|A1|) . 

Moreover, e« = \Jm 2 + p 2 (pi )2 = |P ± p) and P° = E P = v ' M 2 + P 2 are on-shell energies, 
and nucleon and deuteron masses are m and M. This propagator differs from the Dirac 
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two-body propagator of Refs. [[F] only by use of the instant constraint in the frame where 
the two-body system moves with momentum P, instead of in the rest frame. 

As discussed in Ref. ||, it is possible to use the constraint P ■ p = to develop a 
covariant formalism corresponding to instant interactions in the rest frame. If this is done, 
solutions of the equation can be boosted kinematically on a 3D surface embedded in the 4D 
space and defined by the constraint p ■ P — 0. However, this constraint is not compatible 
with the momentum transfer q in interactions. If in the initial state the constraint P ■ 
p = is satisfied, the corresponding constraint for the final state, namely (P + q) ■ (p ± 
|g) = 0, depending on which particle absorbs the momentum transfer, cannot be satisfied. 
Absorption of the momentum transfer requires that the wave function be known off the 
3D surface defined by the constraint. In contrast, a compatible formalism can be obtained 
by adopting the constraint p° = in the Breit frame in place of the covariant one. The 
Breit frame corresponds to initial momentum Pj = (Pp,0,0, — |g) and final momentum 



Pj = (Pp,0,0,|g), where Ep = ^M 2 + q 2 /A. It has the special property that q° = 
for elastic interactions and thus is consistent with conservation of the four-momentum and 
three-dimensional wave functions. 

The inverse propagator is easily obtained from the projection property and it is, 

<7o J (p; P) = 7i°7 2 ° [(Ep/2 - p 1 e l )p 2 + (E P /2 - p 2 ^)pi] , (3) 

with f>i = h(pi)/ei (thus /3jU^(±Pj) = ±uf(±Pj)). The normalization condition for the 
two-body wave function is, 



2E P J (2n 



■d 3 pd 3 k- 

w; P) 



7 o 7 oP 1 +P2 (27r) s 5(p _ k) _ y^(p,k,p/ 



V>(k;P). (4) 



III. BOOST OF THE QUASIPOTENTIAL 

In quasipotential approaches, the quasipotential kernel is formally related to the Bethe- 
Salpeter kernel by, 

K QP (P) = K BS (P) + iK BS (P)A(P)K QP (P), (5) 
where a four-dimensional integration is implied over relative momentum, p, and 

A(P)=G BS (P)-G$ P (P), (6) 

g bs^ = [Pi ■ 7i + m] [p 2 ■ 7 2 + m] 
\p\ — m? + irf\ [p 2 . — m? + irf\ ' 

and 

G$ P (P) = - igo (p,P)2n5(p ) (8) 

in our approach. The quasipotential V used in Eq. ([]]) corresponds to K® p with initial and 
final momenta restricted to the constraint space, i.e., po = 0. As emphasized in the notation, 
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the quasipotential depends in general on the total four-momentum, P. The quasipotential 
propagator involving p° = corresponds to different constraints for different values of P. 
Only when the constraint is expressed covariantly does it have the same physical meaning 
at all P values. 

In order to obtain the quasipotential corresponding to an instant constraint at different 
four-momenta, one must in general solve Eq. (||) at each value of P. Alternatively, it is possi- 
ble to eliminate the Bethe-Salpeter kernel to arrive at a direct relation of the quasipotential 
corresponding to two different momenta P and P . We do this for a somewhat simplified 
case where the Bethe-Salpeter kernel is assumed not to depend on the total momentum, P, 
and find 

K Q (P) = K QP (P ) + tK QP (P ) [A(P) - A(P )] K QP (P). (9) 

In the present work, we consider the boost of the quasipotential from the rest frame, where 
momentum is Pq = (M, 0), to a frame where the four momentum is P = (Ep, P). 

The chief complication in solving for the quasipotential lies in the implied 4D integration. 
However, the problem is soluble for the case of a separable Bethe-Salpeter kernel of the form, 

K=\ X )k(x\, (10) 

where \x) carries the dependence on relative momentum, p, and k is a 16 x 16 matrix. For 
the discussion of the separable potential case, we omit parts of the quasipotential propagator 
which arise from the treatment of crossed Feynman graphs because these contributions are 
not meaningful for a separable potential. A simpler quasipotential propagator is used which 
is expressible as, 

G QP = 2n5(p°) [ ^G BS , (11) 

J ATI 

which agrees with Eqs. and (H) in the + + and rhospin states, but is zero in the 

H — and — h rhospin states. 

It follows that the quasipotential also takes a separable form 

K<3 p (P) = \x)kQ p (P)(x\. (12) 

where k9 p (P) is also a 16 x 16 matrix in the two-particle Dirac space. It is related to the 
Bethe-Salpeter kernel by the matrix equation, 

k QP (P) = k + kA(P)k QP (P), (13) 

where 

d 4 p 

W) 

is a 16 x 16 matrix. Alternatively, one may relate the quasipotential matrices at two different 
momenta by use of Eq. which leads to the matrix equation, 

kQ p (P) = k QP (P ) + k QP (P )[A(P) - A(P ))k QP (P). (15) 

Either Eq. (|T3| ) or (|15| ) is readily solved, e.g., Eq. ( pT3|) yields 

(P) = {1 - kAiP)}}- 1 ^ (16) 



A(P) = • / 7^IX(P) 2 A(P) (14) 
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A. Analysis for a scalar, separable potential 



In order to gain insight into the nature of the boost of the quasipotential, we have solved 
for k QP (P) for the case of a deuteron bound by a scalar, separable interaction. This means 
that k is a coupling constant, g, times the direct product of unit matrices in each particle's 
Dirac space. Taking symmetries into account, we find 

A(P) = (\P ■ 7! +m)(\p -y + m)A l( P) - ^1^A 2 (P) - ^1^A 3 (P) 

( P ■ 7l P ■ 7 2 p. 7l P-7 2 \ 
- 1 M M + ^-^-J A 4 (P) + 71± ■ 72± A 5 (P), (17) 

where P = (£ P , 0, 0, P), P = (P, 0, 0, P P ) and thus P • P = 0. Vectors with the JL subscript 
have only x- and y-components and they are orthogonal to P and P. The quantities A n (P) 
for n = 1 to 5 in this expression are determined by the following equations, 



K{P)=A B n s {P)-^ p {Pl (18) 



where 



A B n S {P) = i J ^ [x{p)f [Pi -m 2 + ivr'lpl - m 2 + i^O^ (19) 
A " P(P) = i I [X(P)]2 I^ [pl - m2 + - ^ + ^" l0 "' (20) 



and 



Oi = 1, (21) 
0l =f^^Y (24) 



M My' 

5 = £ (25) 

Symmetries in the Bethe-Salpeter case cause Af s = Af s , and Af s = 0. 

It is convenient to take matrix elements of A(P) between plane- wave Dirac spinors 
depending on the total momentum and defined as follows, 

u + (F) = N P ( ff -P V (26) 
\E P + M J 

( <J P \ 



«-(P) = iV P I £ P + M j . (27) 
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These states are eigenfunctions of 7 • P/M = ±1 and the normalization factor is 



N P = 



lEp + M 
2M 



(28) 



where Ep = \fP 2 + M 2 . Negative energy states have negative norm: u (P)u (P) = — 1. 
The matrix form of A(P) is, 



[A(P)] 



where 
A 



/A^' ++ (P) _A CT (P) ^A 4 (P) a lz A A (P) \ 

A CT (P) A ' (P) _a lz A 4 {P) a 2± A 4 (P) 

a 2 ,A 4 (P) cx lz A 4 (P) A^' + (P) _A CT ^P) 

V a l2 A 4 (P) a 2 ,A 4 (P) A CT (P) A + ' + {P)J 



\P) = tif(P)ti+(P)A(P)u+(P)ii+(P) = { l -M + mfA^P) - A 2 (P), 



(29) 



(30) 



A ' ~(P) =«r(PK(P)A(PK(PK(P) = (-Af-mfA^P)- A 2 (P), (31) 



A + -' + ~(P) = <(P)n 2 -(P)A(P)<(P)n 2 -(P) = (^M 2 - m^A^P) - A 2 (P), (32) 



A + ' + (P) = «r(P)« 2 + (P)A(P)«r(P)« 2 + (P) = - m 2 )Ai(P) - A 2 (P), (33) 



and 



Act = -a lz a 2z A 3 (P) + <r u • <r 2± A 5 (P) . (34) 

It follows that the solution of Eq. ( |16|) may be found and thus the full structure of 
the quasipotential displayed: k QP = {1 — (yfA} -1 ^! leads to the matrix equation [U — 
sA][l][^ p ]=s[ll],thus, 



[k e > p \ = (g- 1 [1]-[H] [A] [11] 



(35) 



A subtlety here is that the matrix for the unit operator is [11] = diag{l, 1,-1,-1}. It is 
straightforward to realize k® p numerically by calculating the inverse implied in Eq. (j35|). 

In the case of a separable potential, one may readily solve the Bethe-Salpeter wave 
equation or the quasipotential wave equation, using kernels appropriate to each. For the 
Bethe-Salpeter equation, we find 



( j> BS (P)= l ( X \G BS (P)\x)-k- ( j ) BS , 
and for the quasipotential equation, 



(36) 
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2 -2 

P (fm ) 



FIG. 1. Renormalization factor of scalar quasipotential. Solid line shows A for invariant 
deuteron mass and dotted line shows prediction of Eq. (|38|). Dashed line shows result for a scalar, 
meson-exchange interaction. 



QP (P) = i( X \G QP {P)\x) ■ k QP {P) ■ <P QP . (37) 

These are equivalent if k® p is obtained from Eq. fll~3|). The quantity i(x\G B (P)\x) has the 
same form as A(P) using the A PS in place of A n (P). Similarly, i(x\G® p (P)\x) has the same 
form using A® p in place of A n (P). 

Numerical calculations for the scalar separable potential have been performed based on 
using xip) — (1 +P 2 /f J ' 2 ) 2 with /z = 200 MeV for the separable potential. The coupling 
constant g is determined by the condition that a bound state exists for M equal to the 
deuteron mass in the Bethe-Salpeter equation. Solutions corresponding to differing values 
of P are obtained from both equations and they demonstrate that the deuteron mass is 
invariant when the boosted quasipotential is used in Eq. (]37|). 

Owing to dominance of the + + states in the case of weak binding, an approximate 
characterization of the boost of the quasipotential is possible. A significant part of the effect 
is to renormalize the positive-energy matrix element in comparison to its value in the rest 
frame. The ratio of matrix elements defined as in Eq. (|30|) and calculated in the rest and 
moving frames is, 

_ \Tr {k<3 p (P )} ++ ' ++ 
KP) ~ \Tr [kQ p {P)] ++ < ++ ' (38) 

where the trace is over spins. This ratio may be used to approximate the boost as a similar 
renormalization of all matrix elements, 

k QP (P) « k QP (P )/X(P). (39) 



A "fit" renormalization parameter can be determined as a function of P by using Eq. fl39|) 
and solving Eq. ( p7[ ) for the value of A(P) such that the mass is invariant. Figure [1| shows 
the variation of A(P) with momentum for this fit case. For comparison, we show the "ratio" 
prediction based on Eq. (|3lf), which is in reasonable agreement with the "fit" value. Also 
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shown in Figure p] is A(P) fit to yield an invariant deuteron mass with the propagator of 
the next section, but using only the (scalar) a meson. The fit renormalization parameter 
A(P) decreases for increasing P 2 in a qualitatively similar fashion for scalar potentials of 
either the separable or one-boson-exchange type. Equation (|35| ) would remain diagonal and 
proportional to diag{l, 1, —1, —1} if the Dirac structure remained scalar in the boost. As 
|P| increases, the matrix can be seen to change its structure. Owing to the dominance 
of the ++, ++ matrix elements, a significant part of the effect is a renormalization of the 
interaction. 

If one has knowledge of the rest frame interaction /c^ p (P ), Eq. flT5| ) may be used to 
determine k9 p (P). This case is interesting because the NN interaction is usually regarded 
as known in the rest frame and the problem is to boost it to other frames. A more ac- 
curate approximation than a simple renormalization of all matrix elements is to expand 
perturbatively as follows, 

k Q (P) « k QP (P ) + k QP (P ) [A(P) - A(P )] k QP (P ). (40) 

Keeping the second order term in Eq. (20) produces a form for k® p (P) that eliminates most 
of the momentum dependence of the mass. If a renormalization factor A(P) is used with the 
approximation of Eq. (f20|), then a value A = 0.984 at P 2 = 50 fm~ 2 is required to keep the 
mass invariant, as compared with a value A = 0.882 when k® p (P) ~ k® p (Po)/ \ is used and 
A = 1 when the exact k® p (P) is used. 



B. Effective boost approximation 

For the one-boson-exchange potential (in the rest frame) we use the Bonn potential 
(Bonn B, energy-independent, Thompson propagator) fITUl2| . This potential includes scalar, 
pseudo- vector, and vector meson exchanges (a, 5, r], it, lj, p) and is detailed in Appendix [FJ. 
When projected onto positive-energy plane-waves {p\p 2 = ++) in the center-of-mass frame, 
the Dirac two-body go of Eq. (||) reduces to the Thompson propagator and V reduces to the 
Bonn potential. When negative-energies are included, there arise couplings in the quasipo- 
tential for which initial state rho spins of both particles are opposite to the final state rho 
spins. For example, V ++ ' . An analysis of Feynman diagrams where such couplings arise 
shows that they generally are suppressed strongly in comparison with other couplings owing 
to the necessity of large time-like momentum transfer of order 2m. For all other couplings, 
the instant constraint provides a reasonable starting approximation. To accommodate this 
fact, we eliminate the suppressed couplings by setting them to zero, i.e., [j®j®V] Pl ' p 2' pl ' P2 — > 
when p[ = —pi and p' 2 = —p 2 . 

To obtain the correct deuteron binding energy when negative energy sectors are included, 
we modify the Bonn B potential by increasing the scalar attraction about 6%, from g 2 /47r = 
8.0769 to 8.5503. (See Table 0.) 

For the meson-exchange interaction, the simple approximation discussed above has been 
used to boost the quasipotential: V(p' — p, P) = V(p' — p)/A(P), where A(P) is fit to 
produce the correct deuteron total energy, E P = (M 2 + P 2 ) 1 / 2 , using the two-body equation 
and propagator of Eqs. (|IJ) and @. This approximation can be shown to be consistent 
with current conservation when used in the analysis of elastic form factors. Note that the 
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FIG. 2. The scaling of the potential, V(p' — p, P) = V(p' — p)/X(P 2 ), that produces constant 
deuteron mass, M = (Ej, - P 2 ) 1 /2 = 2m - 2.22464MeV: full propagator (solid), ++ states only 
(dotted). 

Dirac spinors appropriate to a moving deuteron are treated exactly and that the factor A(P) 
approximates only the additional change of the potential required in an instant formalism. 
Figure ^ shows that the required change of the potential is modest, with A varying linearly vs. 
P 2 over a wide range of values. For the full meson-exchange interaction, the renormalization 
factor A(P 2 ) increases with P 2 , in contrast with the result of Figure [I] for a scalar interaction. 
This is caused by the differing boost factors required for different types of interaction. When 
A(P) = 1 is used, the potential is too attractive and the binding energy of the deuteron 
increases from 2m - M « 2.2MeV at P 2 = to 2m - M » 4.2MeV at P 2 = 50/m~ 2 . 
We have calculated the deuteron form factors based on A(P) = 1 for comparison with those 
based on A(P) determined so that M is invariant. The differences are small. 




IV. SYMMETRIES AND REDUCTION OF THE EQUATION 

The homogeneous equation is symmetric with respect to the operations of spatial reflec- 
tion (V), particle exchange (JI2), and time-reversal (T), i.e., g = Vg^V^ 1 = V^goP^ 1 = 
TgoT^ 1 , and V is assumed to behave similarly. Thus it is possible to have solutions of the 
homogeneous equation with good corresponding parities: 

^(p,P)=^j / (p,P), (41) 
^(p,P)=r? 12 P 12 ^f(p,P), (42) 
^(p,P) = ^(-l) M T^7 M (p,P), (43) 

where the parity, exchange parity, and anti-linear time-reversal operators are, 

V = 7l 7 2 °^, (44) 
TL2 = V ex U, (45) 
T = afallC, (46) 
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where K, is the operator of complex conjugation and /C and V s reverse three-momenta (V s p = 
— p, V S P = — P); II reverses relative four-momentum (Up = —p); and V ex exchanges Dirac 
indices (V ex j% = JiV ex , V ex u{ 1 (piPi) = u^ 1 (p 1 p 1 )V ex , etc). For the deuteron r] = r) 12 = T]t = 
+1. Note that VT\i is equivalent to the H-parity operator of Kubis Jl3| . 



To find wave functions with good exchange parity, we rewrite the homogeneous equation 
using ip(— p, P) = rii2V ex ip(p, P) and find that, 

V>(p'; P) = g (p'; P)ih°2 I 7^7i°7 2 ° [V(p', p; P) + V(p', -p; P)nuPj[ V>(p; P). (47) 

Note that the necessary range of integration of p z is halved. To find wave functions with 
Mj = and good combined VT parity, we will form eight basis functions with +VT parity, 
and eight with —VT parity. 

The Breit frame total angular momentum operator is, J = Ji + J2 = C + S, where 
£ = 1 + L, 1 = r x p, L = RxP, and S = \(<T\ + c 2 ). Because J 2 and J z commute with 
g^ 1 and V, solutions of the homogeneous equation are eigenfunctions of J 2 and J z . Also, 
wave functions with polarization states Mj = ±1 can be obtained from the Mj = state 



using raising and lowering operators, ip Mj±1 = y (J + Mj)(J — Mj + 1) J ±r i\) Mj . However, 
1, L, and S do not separately commute with g^ 1 , and the usual LSJ partial- wave analysis 
is inapplicable. To proceed, we define cylindrical eigenfunctions of J z which form a basis for 
the dependence of the wave function: 

yZ&fi = e KM ^' l -* ) +\*i)\s2), (48) 

where p = (p xy cos(<f)),p xy cos(0),p z ) is the relative three- momentum, J z y Mj = Mj~y M \ 
and Si = ±| are the +z components of spin (parallel to P = P z z). A related set of 
eigenfunctions which we call the PT basis is given by, 

1 



y% j = ^(y£-y?t), m 

where the ± subscripts are shorthand for ±|, ±|. It will be shown in Appendix [A| that these 
eigenfunctions have good VP parity if Mj = 0. Each set of eigenfunctions is orthonormal, 

l^[yai^y a (<p) = 5 <a , (50) 

where a = (si, s 2 , Mj) (or a = (a, Mj) with a G {el, ol, eO, o0}) and a' is similarly defined. 

To form a complete set of two-particle basis functions, the angular eigenfunctions are 
combined with Dirac spinors obeying the hermitian normalization, u p \pp)u p (p'p) = S pp r 
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X ^(p,P) = <(piPl)^ 2 (P2P 2 )^a(0)- 

Either set of sixteen basis functions is an orthonormal set, 

x^(p,p)f x r P2 (p,P) = 5 P> A>A>, 

and the wave functions are expanded in either set as follows, 

^f( p; p)= x: x£ lrf *(p,PW PB (p*.P n ';P)- 



2tt 



(51) 



(52) 



(53) 



pi,P2,a 



Using these plane-wave basis functions, the homogeneous equation (|T| or ^7) can be 
written in component form as 



k xy dk xy k z 
(2vr) 2 



G p y'\v\v xy \ P)vfif' 2 ' pup2 (p z ,p xy ; k z , k xy )^ p ^k z , k xy ; P) 

(54) 



or 



k z >o (2-n) 



^:>^v\v xv -, k z , k xy w^ p \k\ k xy - p) 



(55) 



where the diagonal hermitian propagator is 



2tt 



xi A (P, P)l <?o(p; P)7!°7 2 ° [X P1 ' P2 (P, P)] = G^V,^; P)^, P A>A',a, (56) 



with 



G p ^(p z ,p xy ;P) 
and the hermitian potential operator is given by, 



(Pi +P 2 )^- -Pip 2 ( e i + e 2 ) 



(57) 



V a p ^ Pl ' P2 (p\p xy ;k z ,k xy ) 



x£'" 2 (p,P) 7i7 2 V(p,k;P) X ^^(k,P) (58) 



or 



(27T) 



x 7l °7 2 ° 



X^" 2 (p,P) 



t/(p,k;P) + t/(p,-k;P)7/ 12 7? s X^ llP2 (k,P). (59) 



A detailed partial-wave analysis of the potential is contained in appendix |B[ Equation (|54"|) 
or (|55|) is solved for Mj = at fixed values of total momentum using the Malfliet-Tjon itera- 
tion procedure |]T4[ and numerical integration over p xy and p z (or radial and polar angle com- 
ponents). Wave functions with polarization states Mj = ±1 are obtained from the Mj = 
state by using the raising and lowering operator, i/j Mj±1 = «/( J + Mj){J — Mj + 1) J ±r i/j Mj . 
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pip 2 lS 


0.0025/m" 2 


p2 

12.5/m~ 2 


25/m" 2 




50/m~ 2 


++01 


0.9499 


0.9433 


0.9363 




0.9226 


++21 


4.97x10 


5.59xl0~ 2 


6.18x10" 


-2 


7.22x10 


+-10 


3.86xl0~ 5 


4.00xl0" 5 


4.09x10" 


-5 


4.18xl0" 5 


+-11 


1.46xl0~ 4 


1.19xl0" 4 


9.80x10" 


-5 


6.86xl0~ 5 


01 


2.70xl0~ 7 


3.07xl0~ 7 


3.20x10" 


-7 


2.82xl0~ 7 


21 


2.49xl0~ 6 


1.72xl0~ 6 


1.27x10" 


-6 


8.01 xlO" 7 


P{++) 


0.99962 


0.99921 


0.99815 




0.99475 



TABLE I. Wave function projection probabilities (J = 1 and Mj = 0). See text. 
P(++) = P(++01) + P(++21). 

V. RESULTS FOR WAVE FUNCTIONS 

Our wave functions vary as a function of both the magnitude and polar angle of relative 
momentum. To show how the wave functions change with total momentum, we project 
them onto standard LSJ basis functions and integrate out dependence on the polar angles, 

«k(|p|> p) = / dn p [y lsJ M.M ^ Mj (p, p), (so) 

where the LSJ basis functions are given by, 

ylsj(P)= E (l,S,m l ,m s \l,S,J,M J )Y l>mi (d,<i ) )\S,m s ), (61) 

m s =-l 

where I and S are the relative-orbital and spin angular momenta, j = 1 + S and 1 = r x p. 
These wave function components are shown in Fig. |] for total momentum corresponding 
to P 2 =.0025, 12.5, 25 and 50/m~ 2 . A 'probability' for the wave function projections as 
defined by 

P iT 2 = ^Jp 2 d P \^ P2 (\p\,P)\ 2 (62) 
is shown in Table | based on the wave function normalization of Eq. (|J). 

VI. CONCLUDING REMARKS 

The theory of relativistic bound states is formulated in three dimensions by use of an 
instant reduction of the Bethe-Salpeter formalism. This formalism is applicable to the 
analysis of elastic form factors and it is consistent with current conservation and four- 
momentum conservation. Elastic form factors involve matrix elements which are calculated 
in the Breit frame, thus requiring wave functions for the initial and final states that have 
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FIG. 3. Dependence of Breit-frame wave functions on total momentum. In this figure 
J = 1 and Mj = 0, P » OMeF(solid), P = 698MeF(dotted), P = 987MeV (dashed), and 
P = 1395MeF(long dots). 
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been boosted to momentum P = ±|q. These wave functions are calculated in this paper for 
the case of the deuteron. The form factors for the deuteron will be the subject of another 
paper. 

The main issues addressed in this work are the boost of the interaction that is required 
and the solution of the quasipotential wave equation in frames where the total momentum is 
nonzero. We have shown that in general the boosted interaction is defined as the solution of a 
four- dimensional equation and that equation has been solved for the special case of a scalar, 
separable interaction. The main effect of the boost is to renormalize the dominant matrix 
elements of the interaction as three-momentum varies, although there is also in general a 
change in the Dirac structure of the interaction. An approximation which captures the 
renormalization effect is used for the more complicated one-boson exchange interaction. For 
momenta up to P 2 = 50 fm -2 , the renormalization of the interaction is modest in the case 
of the deuteron, varying linearly with P 2 and amounting to about a 10% reduction of the 
one-boson exchange potential at P 2 = 50 fm~ 2 . A complete solution of the boost of the 
meson-exchange interaction is left as an unsolved problem. However, the separable potential 
analysis suggests that a perturbative expansion of Eq. fll5]) may provide accurate results for 
the boost. 

The solutions for the quasipotential wave functions have been developed and the results 
show that there are modest variations as the deuteron is boosted to momenta up to P 2 =50 
fm -2 , corresponding to q 2 = 200 fm -2 . In a future article, they will be applied to the 
calculation of elastic form factors for the deuteron. 

Support for this work by the U.S. Department of Energy under grants DE-FG02-93ER- 
40762 and DE-AC05-84ER40150 is gratefully acknowledged. 



APPENDIX A: BASIS FUNCTIONS AND OPERATORS 

The hermitian plane-wave spinors for particle i and denoted by (p^k) are 




(Al) 



with N = f^2e~) ^ anc ^ €i = ( m2 k?) 1 ^ 2 - It will be convenient later to use spinors with 
zero momentum which are related to the plane-wave spinors by, 

uf ( P M) = N(l + Pi ^>f (0). (A2) 

Symmetry and Dirac operators acting on the basis functions will produce linear combi- 
nations (independent of 0) of the basis functions. Some useful examples of this follow. The 
action of the operators on the Dirac and spin parts of the basis functions will be considered 
separately. Consider first the symmetry operators of Eqs. (041) - (|46l), 
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Puf(piki) = Pi uf{ Pi \n)V s , (A3) 
P 12 «f ( P X) = uf{ Pi k,)V l2 (j ^ i), (A4) 
Tu?(pX)=<( Pt k l )T, (A5) 



and 



v s y^M = (-i) Mj - si - S2 y^M> (A6) 
vuy&w = (-i) Mj - si - S2 y^ sl (<P), (at) 
r^fo) = (-i) Mj ^_ S2 (0). (as) 

Similar relations for the 3^* fj (0) of Eq. (fj9|) (a G {el, ol, eO, oO}) are easily derived from the 
above equations. In particular, for the basis XMj P «(k, P) = u^ 1 (^iki)^ 2 (p 2 k 2 )y^ Ij (0), 

OTv Pi,p 2 /u p n / +PiP2Xa'-M J (k,P) a = ol,eO 

(A9) 

Thus these basis functions have VT = ±1 if Mj = 0. The condition VT = ±1 is equivalent 
to Parity = V = ±(— 1) J+1 for the standard LSJ basis functions in the center-of-mass frame. 

Consider next Dirac operators such as 7°, 7^ • k i: ji ■ 72, etc., acting on the simple basis 
functions given by x^'g(k, P) = u^ 1 (0)u^ (0)y^f' L (</>)■ Their effect on the zero-momentum 
spinors is particularly simple, 

7 l V l (0) = P,< 1 (0), (A10) 
7?tif(0)=ur*(0). (All) 

Moreover, spin operators acting on the basis functions produce linear combinations as fol- 
lows, 

<xi ■ <T 2 y™ S2 = is lS2 y™ S2 + 25 Sl ,. S2 y M Si ^ (A12) 

<n • = 2 Sl kiy™ S2 + kf y y^ S2 , (Ais) 
<x 2 ■ k 2 y™ S2 = 2s 2 k* 2 y™ S2 + k* 2 y y™- S2 , (au) 

where ki V = k xy = —k 2 y . For the analysis of potential operators, it is useful to note, 

*i ■ nyZsM = zsitiK„{<i>k) + e 2si ^ - ki v )y M 81 ,sM, (AIS) 
°2 ■ ^y^sM = 2s 2 qiyM S2 (<f> k ) + ( P fe 2 -^ - k* 2 v)yM_ S2 (<f> k ), (Aie) 

where q; = pi - k { and <p d = <p p - (j) k . 

Similar relations can be derived for the same Dirac operators acting on the sixteen 
Ui 1 (0)u 2 2 (0)3^q /j (0) basis functions. This can be accomplished by using the above equations 
and defining a basis transform between the two bases. Note that some of the above operators, 
such as 7I, do not commute with VT. 



16 



APPENDIX B: PARTIAL- WAVE POTENTIAL ANALYSIS 



In this appendix we sketch the partial wave analysis of the one-boson-exchange potential. 
This is most readily accomplished using the Xsl% = u i 1 W u 2 > (0)3^fi 2 (</>) basis functions. 
The potential in this basis can easily be transformed to plane-wave or PT bases. The 
operator form of the potential is generated from the Feynman rules for meson propagators, 



A(q) 



and meson-nucleon vertices, 



q 2 — fi 2 + in' 



-9a,P + 



(Bl) 
(B2) 



A*(5) 
Ap(«) 



2m 



(B3) 
(B4) 

(B5) 



where the exchanged four momentum is q = p\ — k\ = p — k = k 2 — p 2 . The exchange 
of a single meson is given by A(q)A(q)A(—q) (or A a (g)A Qi/ g(g)A /3 (— q) for vector mesons) 
with a factor of T\ ■ t 2 added to the exchange of isovector mesons and a coupling constant 
g and form factor F(q 2 ) attached to each meson-nucleon vertex. Note that for the isoscalar 
deuteron, n • r 2 — > —3. The Bonn model form factors are, 



F t (q 2 



A? 



Af - g 2 

Thus the full one-boson-exchange potential is, 

?(p,A0=E*^(?)iW)%), 



(B6) 



(B7) 



where i is summed over the six mesons of Table [H| 



n -t 2 if i 

i z-/ i 



isovector meson 
isoscalar meson 



and 



(B8) 



t>5(?)=A s (g)A 5 (- g ), 
Vp( ? ) = -Ap(g)A P (-g), 



Vy(?) = A a v (q) 



Using Eqs. (|ATo|) to flATop it can be seen that 



-<i)- 



Xs",s"Mj 



;k,p)[v^(^,^,F,^)] 



s' 1 ',S2,si,S2,A^j;n' : 



(B9) 
(BIO) 

(Bll) 



(B12) 
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Bonn B potential parameters 
(energy independent, Thompson propagator) 





7T 


V 


P 


UJ 


5 


a 




PV-IV 


PV-IS 


V-IV 


V-IS 


S-IV 


S-IS 


n 


.13803 


.5488 


.769 


.7826 


.983 


.550 


A 


1.2 


1.5 


1.3 


1.5 


1.5 


2.0 


9 2 /4vr 


14.6 


5.0 


.95 


20.0 


3.1155 


8.0769 


/ 






6.1 









g 2 /47r with ne; 


native enerj 


ry sectors 








8.5503 



TABLE II. Meson Parameters of the Bonn B model [12|, and modification to g 2 /47r when neg- 
ative energy sectors are included. M and A in GeV. Nucleon mass is m = .938926Gey, deuteron 
binding energy is 2.224644MeV. Vector propagator is approximated by A a ^(q) — > —g a ^A{q). 



where i = S, P,V, n = 0, ±1, ±2 and <pd = 4>p — 4>k with implicit sum over double prime 
variables. The exact form of [V^] is easily derived using a symbolic manipulation program 
such as Mathematica. For Vy, note cos(0d) = (e l ™ — e~ l ^ d )/2. Next, the scalar part of the 
potential, which is a function of <pd but not of <f> p or alone, can be integrated, 



(27T) 



where I = Mi 



Xs' v s' 2> M'j [ 



Are 



2 And 



M ( 



rn—l 



<Vi ,p'/ <V 2 ,pg ^ ^4 ,4' ^,Mj JfAF) 



s' 2 and 



jn 

FAF 



f ^A(q)F 2 (q 2 )e in ^. 



Combining the necessary factors, the partial-wave potential is, 



P2 

S2,Mj;n' 



(B13) 



(B14) 



(B15) 



This partial-wave potential based on zero-momentum spinors can be transformed to a plane- 
wave basis using the basis transform implied by Eqs. (|A2| ), ( [All] ), and ( |A13| - [AT4] ), 

Finally, the scalar <ftd integral of Eq. QB14Q must be evaluated. This can be accomplished 
analytically using z = and contour integration around the unit circle in the complex 
plane. We first express AF 2 as a sum of simple denominators times <pd independent coef- 
ficients using a partial fraction expansion. The integral of a simple meson denominator is 
given by, 



d(f) a 



in<j> d 



2n a^-bcos((p d ) J 2n 



-A(q)e m 



4>d 



(B16) 



where a M = $ - (q(<j> d = 0)) 2 , b = 2p xy k xy (thus -A x (g) = /i 2 - q 2 = a u - 6cos(0 d )), 
and a u > b, z u = (a — S)/b, S u = \Ja 2 — b 2 . The denominators of the Bonn form factors, 
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Eq. ( |B6|) , take the same form with // replaced by A; thus I a = I^fi — > A). Using the partial 
fraction expansion for Bonn form factors and the simple integrals above, 

Ifaf = -i; + 11- (A 2 - (BIT) 

This integral can be numerically checked using the following simple, well known, technique. 
We can write, 

A(g)FV) = A(g) + (*V) - F 2 {q 2 )) A(g), (B18) 

for any form factors with F(fj, 2 ) = 1. The first term is integrated using Eq. ( |B16| ), while the 
second term is evaluated numerically. Note the second term is non-singular even at q 2 = /i 2 , 
which can not occur with instant constraints, but can occur if different constraints are used 
on the left and right of V. 
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